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摘 要 随 着 空间 数据 的 大 量 增长 , 对 数据 可 视 化 和 数据 存 取 效 率 提出 了 更 高 要 求 , 迫切 需要 对 数据 进行 有 效 的 
组 织 和 管理 . 对 庞大 的 日 地 空间 , 采用 SDOG-R 方法 将 日 地 空间 剖 分 为 不 同 分 辩 率 等 级 的 格 网 ,并 针对 该 网 格 提 
出 相应 的 编码 方案 .以 太阳 风 模 型 数据 为 例 , 给 出 了 具体 的 组 织 实例 , 经 实验 验证 , 该 剖 分 模型 不 仅 解决 了 球 心 处 
网 格 过 密 问 题 , 还 满足 了 径 向 分 辨 率 大 于 经 纬 球面 分 辨 率 的 需求 ， 基 于 三 维 立 体 剖 分 的 太阳 风 LOD 空间 数据 模 
型 , 不 但 能 提供 多 分 辩 率 数据 , 而 且 显 著 提 高 了 大 规模 数据 检索 和 存 取 速度 , 有 效 地 支持 海量 空间 数据 的 组 织 管理 . 
关键 词 空间 痢 分 , SDOG-R, 编码 , 多 分 辨 率 
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Abstract With the ever-increasing of space data, a reasonable data-arrangement is strongly re- 
duired for high efficiency in data accessing and visualization. In this paper, basing on the character- 
istics of solar-wind data, a SDOG-R, grid model is employed, which is radial independent division in 
the adaptive SDOG grid. For huge Sun-Earth space, it is divided into multi-resolution grids. The 
original data, which is output from a SIP-CESE solar wind model, are re-sampled, and put into 
the subdivided grids. The grid model is applied to both regular and irregular sampling data, and 
fully maintains the characteristics of original data to the best and encodes grids with an improved 
CDZ curve. Specific examples of data, organization are given. The experiments prove that the grid 
model not only resolves the problem that grids are too dense at two poles and the spherical center, 
but also meets the need for higher resolution in radius than in latitude and longitude. Besides, 3D 
LOD spatial data model can not only provide multi-resolution data, but also significantly improve 


the large-scale mass data retrieval and access efficiency, and can support the organization and man- 


2014-07-26 收 到 原稿 , 2015-05-17 收 到 修 定稿 
E-mail: hyssmilyQ163.com 


胡 雅 斯 等 : 基于 日 心 坐标 系 的 三 维 立 体 剖 分 模型 及 编码 


agement of massive spatial data effectively. 
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0 引言 


中 国 针对 太阳 、 地 球 、 月 球 的 空间 探测 开展 了 
一 系列 项 目 研究 ,例如 夸 父 计划 、 地 球 空间 双星 探 
测 计划 、 高 分 状 率 对 地 观测 系统 工程 、 探 月 工程 等 . 
这 些 探测 项 目 已 经 获取 或 即将 获取 到 大 量 空 间 科学 
探测 数据 , 例如 磁 层 、 辐 射 带 和 电离 层 数据 等 . 在 科 
学 研究 中 除 探 测 数 据 外 , 还 需要 提供 一 些 模型 输出 数 
据 , 例如 太阳 风 数 据 等 . 一 方面 , 空间 信息 量 的 不 断 增 
长 并 不 意味 着 有 效 的 空间 信息 必然 增长 , 无 序 的 空间 
信息 资源 有 可 能 加 剧 信息 增长 与 有 效应 用 之 间 的 矛 
盾 上 3; 另 一 方面 , 如 何 从 当前 大 量 复杂 的 数据 中 找 
到 需要 的 数据 是 迫切 需要 解决 的 问题 加 

地 理 信息 系统 的 发 展 , 特别 是 数字 地 球技 术 的 发 
展 , 使 地 球 科学 研究 具有 强 有 力 的 数据 组 织 和 管理 手 
段 . 面 对 快速 发 展 的 空间 科学 , 对 数据 进行 统一 组 织 
管理 , 提高 数据 的 可 视 化 效果 和 数据 存 取 效率 , 具有 
重要 理论 研究 意义 和 广泛 的 应 用 前 景 Bl. 

在 庞大 的 太阳 系 空间 , 为 了 研究 不 同类 型 空间 数 
据 的 变化 情况 , 需要 构建 三 维 立 体 训 分 格 网 模型 并 且 
对 网 格 进行 编码 , 进而 将 数据 划分 成 不 同 的 分 辩 率 等 
级 , 实现 对 庞大 空间 数据 的 组 织 , 提高 存 取 效率 加 ,为 
空间 数据 的 组 织 、 管 理 、 分 析 提供 技术 支撑 斑 2 中 

由 于 目前 中 国 尚 没有 实际 探测 的 太阳 风 数 据 , 本 
文 将 采用 中 国 空间 天 气 学 国家 重点 实验 室 研发 的 三 
维 太阳 风 模 型 输出 数据 作为 测试 数据 15- 坷 , 对 日 
地 至 间 进 行 立体 放 分 , 并 提出 了 相应 的 编码 方案 . 


1 齐 分 模型 的 构建 


目前 , 立体 痢 分 网 格 模型 研究 大 都 基于 地 球 空间 
的 剖 分 , 为 了 处 理 大 量 的 地 学 数据 , 需要 将 整个 地 球 
实体 进行 格 网 划分 , 格 网 模型 有 多 种 , 包括 Yin-Yang 
格 网 9、Stemmer 立体 格 网 PCI、Ballard 格 网 分 层 
结构 2、Stadler 立体 格 网 [3 等 . 正 多 面体 立体 剖 
分 模型 有 些 是 基于 四 面体 和 六 面体 的 , 其 运算 量 大 ， 
通用 性 不 强 . 为 解决 这 些 问 题 , 不 断 有 新 的 格 网 划分 
模型 提出 , 包括 QuaPAP3 、 球 体 退 化 八 又 树 格 网 S- 


DOG (Sphere Degenerated-Octree Grid)24-25、 适 应 
性 SDOG 格 网 Bal. 

适应 性 SDOG-RE9 格 网 模型 是 在 SDOG 格 网 
的 基础 上 , 沿 径 向 递归 进行 二 分 细 分 , 直至 格 网 的 径 
向 分 辩 率 满足 要 求 为 止 ， 该 模型 既 继承 了 球面 退化 
四 又 树 的 优点 , 又 解决 了 球 心 处 格 网 过 密 问题 , 由 于 
在 半径 上 可 继续 进行 独立 剖 分 并 能 满足 径 向 分 辨 率 
大 于 球面 分 辨 率 的 需求 , 因此 本 文选 用 SDOG-R 格 
网 模型 实现 空间 数据 的 三 维 立 体 剖 分 . 
1.1 SDOG-R 格 网 模型 

SDOG-R 格 网 的 具体 放 分 步骤 26 如 下 . 

步骤 1 将 整个 球体 在 经 度 上 平分 2 次 , 纬度 上 
平分 1 次 , 形成 8 个 八 分 体 , 如 图 1(a) 所 示 ， 

步骤 2 八 分 体 第 1 次 剖 分 . 首先 , 等 分 径 线 ， 
以 半径 为 原 半 径 一 半 的 球面 切割 八 分 体 , 将 八 分 体 分 
为 内 网 格 和 外 网 格 , 如 图 1(b) 所 示 ; 然后 , 等 分 经 线 ， 
取 外 网 格外 表面 左右 两 经 线 的 中 点 , 用 纬 线 连接 两 中 
点 , 再 以 新 产生 的 纬 线 与 球 心 所 构成 的 圆锥 面 切割 外 
格 网 , 生成 外 上 网 格 和 外 下 网 格 , 内 格 网 保持 不 变 , 如 
图 1(c) 所 示 ; 最 后 , 等 分 纬 线 , 取 外 下 网 格 的 上 下 纬 
线 中 点 ， 以 两 中 点 及 球 心 确定 的 平面 切割 外 下 网 格 ， 
生成 外 下 无 网 格 和 外 下 右 网 格 , 如 图 1(d) 所 示 . 

步骤 3 八 分 体 首次 剖 分 将 产生 4 个子 网 格 , 根 
据 网 格 不 同 的 退化 特点 , 将 内 单元 格 命名 为 球面 退化 
格 网 (SG), 外 上 单元 格 命名 为 纬 线 退 化 格 网 (LG), 外 
下 左 单元 格 和 外 下 右 单元 格 命名 为 正常 格 网 (NG). 
SDOG 格 网 全 部 由 这 三 种 基本 格 网 构成 . 经 过 第 1 
次 训 分 生成 一 个 SG, 一 个 LG, 二 个 NG. 

步骤 4 在 第 1 层 SDOG 格 网 基础 上 , 进行 第 2 
次 剖 分 , 对 SG, LG, NG 三 种 格 网 分 别 剖 分 , 其 中 SG 
按照 步骤 3 的 规则 痢 分 . 纬 线 退化 格 网 LG 要 经 过 退 
化 , 纬 线 退 化 格 网 LG 的 内 层 格 网 和 外 层 格 网 均 痢 分 
成 3 块 , 则 1 个 格 网 分 成 6 块 ; 正常 格 网 NG 则 是 八 
又 树 剖 分 , 且 1 个 格 网 分 成 8 块 . 

步骤 5 重复 上 述 步 又 2~4, 生成 不 同 分 辩 率 等 
级 的 SDOG 格 网 ( 见 图 2~4). 

步骤 6 在 每 一 层 的 SDOG 格 网 基础 上 , 在 径 
向 继续 细 分 , 得 到 对 应 剖 分 层次 的 SDOG-R 格 网 , 细 
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分 次 数 视 需求 而 定 。 图 5 所 示 为 在 3 次 剖 分 的 
SDOG 格 网 基础 上 , 对 某 一 个 SDOG 网 格 在 径 向 3 
次 细 分 的 结果 , 一 个 SDOG 网 格 则 变 成 8 个 子 网 格 ， 
SDOG-R 格 网 生成 . 

通过 径 向 截面 结合 球面 剖 分 可 以 分 析 SDOG 格 
网 剖 分 过 程 , 进而 给 出 格 网 数量 计算 方法 . 

图 2(a) 给 出 的 是 八 分 体 径 向 剖 分 结果 从 
图 2(a) 可 以 看 出 , 第 1 次 剂 分 时 径 向 1 分 为 2, 从 


(a) eighth sphere 


NG 


(c) subdivide in longitude 
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内 到 外 分 别 为 第 0 层 和 第 1 层 , 图 2(b) 给 出 的 是 外 
侧 第 1 层 对 应 球面 的 网 格 训 分 结果 . 图 2(a) 网 格 上 
的 数字 为 对 应 球面 的 网 格 个 数 , 上 网 格 1 个 则 标 1， 
下 网 格 有 左右 2 个 则 标 2, 第 1 次 放 分 一 共有 4 个 网 
格 . 

如 图 3 所 示 , 第 2 次 齐 分 时 径 向 1 分 为 4. 可 以 
看 出 , 第 2 次 剖 分 其 实 是 第 1 次 剖 分 的 结果 全 压缩 在 
八 分 体内 侧 一 半 , 而 外 侧 一 半 则 是 第 1 次 第 2 层 继 


inner grid 
/ SG 


(b) subdivide in radius 


outer under right grid 
NG 


(d) subdivide in latitude 


图 1 SDOG 格 网 第 1 次 剖 分 
Fig.1 ‘The first subdivision for SDOG 


(a) 


upper grid 
1 


under grid 


number in radius 0 1 


图 2 第 1 次 训 分 SDOG 格 网 数量 
Grids’ number diagram for the first SDOG subdivision 


Fig.2 


(b) 


upper grid 


under left | underright\ 2 


grid 


grids in 1st sphere after the 1st subdivision 
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number in radius 0 1 2 3 


图 3 第 2 次 训 分 SDOG 格 网 数量 


A 


Fig.3 Grids’ number diagram for the second SDOG subdivision 


图 4 第 3 次 训 分 SDOG 格 网 数量 


Fig.4 Grids’ number diagram for the third SDOG subdivision 


续 训 分 的 结果 . 依次 类 推 , 第 3 次 训 分 结果 如 图 4 所 
不 : 

从 图 4 可 以 看 出 , 以 半径 为 一 半 的 位 置 为 界 , 半 
径 由 外 往 里 , 网 格 相 对 密度 是 递增 的 , 从 半径 1/2 处 
到 半径 1/4 处 这 一 段 , 相对 网 格 密度 也 是 递增 的 , 但 
是 这 段 密度 只 有 外 侧 1/2 的 一 半 . 所 以 网 格 密度 在 
半径 1/2 处 附近 达到 最 高 , 其 次 是 1/4 处 . 表明 当 齐 
分 层次 为 n 时 , 径 向 齐 分 得 到 的 球面 层 数 是 2”. 

从 图 4 还 可 以 发 现 , 除去 北极 项 的 1 个 网 格 比 
较 特 殊 外 , 其 他 网 格 都 与 2 有 关 , 其 具有 一 定 规律 性 ， 
第 n 次 剖 分 , 外 侧 半球 被 分 成 2"” 层 球面 , 最 外 层 


的 球面 网 格 个 数 N 计算 如 下 : 
n=1 时 , N=2=1x2=20x2!; (1) 
n=2 时 N=2+2x4= 
20 x 21 + 21 x 22; (2) 
n=3 时 N=242x444x8= 


90 22 让 二 
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(b) 
1 
2 
4 
4 
grids in 2nd and 3rd after 2nd 
sphere subdivision 
grids in Sth 6th 7th 8th 
sphere after 3rd subdivision 
22 3 (3) 


n= 二 4 时 ,N=2 +H2x4++4x8+ 二 8x16= 
20 x 21 +21 x 22 + 22 x 23+4 


23 x 24, (4) 
n =i 时 , N = 2° x 21+421 x 22+ 

22 x 2 十 十 21 文 2 三 

21 十 23 十 25.. .十 22 一 ! 一 

办 

5( 人 一 1 (5) 
第 ;次 训 分 后 , 外 侧 半球 网 格 总 数 为 

2 广 Ay —1). (6) 


如 果 经 过 ?次 放 分 , 则 不 包含 1 的 体 剖 分 格 网 数 为 
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1 Nn n 


420" -2 (>) 本 


= Sl 


se -D)—2(2"—1)|= 
ee —7x2"t +6). (7) 


再 加 上 所 有 未 考虑 的 1, 由 于 每 一 层 顶 上 都 有 一 个 1， 
经 过 呈 次 放 分 , 分 成 2 层 , 则 有 2” 个 1. 因此 , 八 分 
体 所 有 网 格 个 数 为 


N=N+2"= 
下 
re —7x2"t!l+6)+2". (8) 


对 于 普通 八 又 树 剖 分 得 到 的 网 格 个 数 N = 23” (n > 
0), 可 以 发 现 格 网 数量 以 23 速度 增长 . 
表 1 给 出 了 经 5 次 剖 分 的 网 格 数量 统计 . 从 表 1 
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可 以 看 出 普通 球体 八 又 树 网 格 与 SDOG 网 格 的 比较 
结果 当 痢 分 层次 大 于 3 的 时 候 , 数据 压缩 比 就 稳 
定 在 61%. 这 样 既 节 省 了 大 量 的 数据 存储 空间 和 格 
网 编码 数量 , 又 有 效 平衡 了 两 极 格 网 过 密 的 缺点 , 使 
得 格 网 大 小 趋 于 均匀 化 ， 然而， 当 径 向 分 辨 率 大 于 
球 向 分 辨 率 时 , SDOG 格 网 已 不 能 满足 需求 , 为 了 最 
大 程度 保留 径 向 数据 的 完整 性 , 解决 两 个 维度 上 分 辨 
率 不 一 致 的 问题 , 可 以 在 SDOG 格 网 的 基础 上 继续 
在 径 向 独立 细 分 , 使 得 径 向 的 剖 分 次 数 大 于 球面 的 
剖 分 次 数 , 以 满足 径 向 分 辨 率 大 于 球面 分 辨 率 需 求 . 
因此 , 提出 径 向 细 分 的 球体 退化 八 义 树 格 网 (Sphere 
Degenerated Octree Grid-Radius), 即 SDOG-R 格 网 . 

综 上 所 述 , 每 一 个 SDOG-R 格 网 都 是 在 对 应 剖 
分 层次 的 SDOG 格 网 基础 上 生成 的 , 其 中 径 向 细 分 
的 次 数 要 根据 实际 需求 进行 相应 调整 . 特别 是 当 径 
向 细 分 次 数 为 0 时 , 即 是 SDOG 格 网 . SDOG 格 网 


图 5 第 1 层 SDOG 格 网 径 向 3 次 细 分 构架 
Fig.5 The third subdivision in radius for the first SDOG 


表 1 八 分 体 剖 分 网 格 数量 比较 表 


Table 1 Comparison table for grids number after subdivision of eighth sphere 
痢 分 “普通 球体 八 又 树 。SDOG 格 网 数量 
层次 ” 格 网 个 数 格 网 数 降低 比率 /(%) 
Tr “下 二 站 1 十 2 十 1=4 50.00 
2 25=64 [1+2)x2+2x8]+4=26 59.38 
3 29=512 [1 x6+(2+4+4) x8]x2+26=198 61.33 
4 2!'?=4096 [1x6+(2+4x2+8x4)x8]x4+198=1566 61.77 
5 2 =32768 [1 x6+(2+4x2+8x4+16x8)x8]x8+1566= 12494 61.87 
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主要 是 用 较 少 的 网 格 数 解决 质心 格 网 过 密 的 问题 
而 SDOG-R 主要 是 解决 径 向 和 球面 两 个 维度 分 辩 率 
不 同步 的 问题 . 
1.2 ”SDOG-R 格 网 编码 

编码 是 将 每 个 谢 分 后 的 子 网 格 赋予 唯 一 的 标识 ， 
使 得 网 格 中 的 坐标 、 属 性 数据 和 编码 形成 一 一 对 应 
的 关联 关系 . 为 了 进行 高 效 的 数据 索引 , 格 网 编码 应 
遵循 相 邻 格 网 的 编码 尽 可 能 连续 ， 连 续 的 编码 所 对 
应 的 格 网 尽 可 能 相 邻 的 原则 ， 下 面 以 第 1 层 放 分 
SDOG-R 为 例 , 详 述 其 编码 规则 24271. 

(1) 用 0~7 的 八进制 数 对 八 分 体 进行 编码 , 以 确 
定 八 分 体 所 在 的 象限 ， 按 照 逆 时 针 方向 将 北半球 依 
次 编号 为 0, 1, 2, 3, 将 南半球 依次 编号 为 4, 5, 6, 7， 
如 图 6 所 示 . 

(2) 参照 图 1, 八 分 体 的 首次 前 分 即 沿 径 向 1 分 
为 2, 将 内 层 记 为 0, 外 层 记 为 1; 每 层 都 分 成 4 个 网 
格 , 0 层 4 个 网 格 合并 为 1 个 , 用 0 表示 , 即 SG 网 
格 ; 1 层 上 面 的 两 个 网 格 合 为 一 个 , 用 4 表示, 即 LG 
网 格 ; 下 面 两 个 网 格 即 左下 网 格 , 右 下 网 格 , 分 别 编 
码 为 6, 7, 即 NG 网 格 . 

(3) 对 于 八 分 体 第 2 次 训 分 , 顶 上 三 楼 柱 编码 与 
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第 1 次 剖 分 单元 的 编码 方法 相同 , 一 分 为 六 , 新 网 格 
编码 在 之 前 的 基础 上 分 别 加 上 0, 2, 3, 4, 6, 7; 子 六 面 
体 一 分 为 八 , 新 编码 分 别 在 之 前 的 编码 基础 上 加 上 0， 
1, 2, 3, 4, 5, 6, 7, 依次 类 推 ; 剖 分 层次 每 增加 1, 编码 
相应 增加 一 位 . 

(4) 分 割 线 后 面 编 码 标识 的 是 径 向 细 分 后 的 层 
号 , 用 二 进 制 码 表 示 .， 在 径 向 上 ， 每 个 SDOG 格 
网 都 要 经 过 径 向 3 次 痢 分 , 即 一 分 为 八 , 则 任意 一 
个 SDOG-R 格 网 细 分 后 , 在 其 本 身 编 码 后 加 上 径 向 
细 分 标识 位 , 径 向 细 分 标识 位 确定 SDOG-R 格 网 的 
具体 位 置 . 由 图 7 可 以 看 出 , 在 第 一 次 SDOG 网 格 的 
基础 上 , 继续 进行 径 向 1 次 细 分 , 分 割 线 后 面 的 编码 


latitude 
90° 


longitude 
十 180° 


图 6 八 分 体 编码 
Fig.6 Coding for eighth sphere 
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Sphere 1 


Sphere 0 


| We SOG 
六 Ee 
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in radius 
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图 7 第 1 层 


层 SDOG 径 向 1 次 细 分 编码 


Fig.7 Encoding diagram for the first subdivision in radius for the first SDOG 
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标识 是 径 向 细 分 后 的 层 号 , 用 二 进 制 码 表示 . 径 向 3 
次 细 分 后 的 结果 如 图 8 所 示 . 综 上 所 述 , 要 唯一 标 
识 一 个 SDOG-R 网 格 , 编码 包含 两 部 分 , 即 分 割 线 前 
的 SDOG 格 网 编码 和 分 割 线 后 的 径 向 层 号 . 编码 可 
以 表示 为 

qngqn+1—P. (9) 


式 中 , qi 为 八 分 体 象限 标识 位 , q; (i > 1) 为 每 一 层 八 


Morton = gqg192 …， 
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分 体 放 分 每 个 网 格 对 应 的 编码 , 符号 “一 ”为 分 隔 符 ， 
Pp 为 以 SDOG 格 网 为 基础 在 径 向 上 细 分 3 次 后 的 径 
向 标识 位 . 

径 向 标识 位 p 与 原 SDOG 编码 用 分 割 线 分 开 ， 
2 用 二 进 制 码 表示 , 1 位 代表 细 分 1 次 , 2 位 代表 细 
分 2 次 , 每 经 过 一 次 细 分 , 编码 增加 一 位 , 剖 分 的 次 数 
越 多 , 生成 的 子 网 格 越 小 . 依 此 类 推 , 第 2 层 SDOG- 
R 格 网 编码 如 图 9 所 示 . 


AO 


图 8 第 1 层 


3rd subdivision 
in radius (SDOG-R) 


层 SDOG 径 向 3 次 细 分 编码 


Fig.8 Encoding diagram for the third subdivision in radius for the first SDOG 
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图 9 第 2 层 
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Fig.9 Encoding diagram for the first subdivision in radius for the second SDOG 


胡 雅 斯 等 : 基于 日 心 坐标 系 的 三 维 立 体 剖 分 模型 及 编码 


通过 网 格 编码 , 使 得 网 格 中 的 数据 与 编码 相互 对 
应 , 为 数据 检索 、 组 织 及 数据 应 用 与 分 析 提 供 技术 文 
持 . 


2 数据 组 织 实例 


以 太阳 风 粒 子 密度 数据 为 例 ， 阐述 如 何 使 
用 SDOG-R 模型 对 空间 数据 进行 组 织 和 管理 ， 该 
数据 是 由 太阳 行星 际 守 恒 元 解 元 (SIP-CESE) 三 维 
太阳 风 模 型 06- 5 输出 的 . 

通过 对 该 太阳 风 数 据 的 组 织 , 通过 可 视 化 展示 可 
以 帮助 科学 而 家 分 析 该 太阳 风 模 型 是 否 满足 预期 , 能 
为 科学 研究 提供 技术 支撑 [74. 

2.1 ”数据 分 析 

所 分 析 的 数据 在 极 坐标 下 采样 ， 三 个 坐标 轴 分 
别 为 半径 、 经 度 和 纬度 ， 其 具有 不 规则 采样 的 特 
点 ， 在 以 太阳 为 中 心 的 球体 采样 空间 中 , 纬度 范围 
为 -90”~ 90° 采样 55 次 , 经度 范围 为 0>~360° 
采样 80 次 , 径 向 采样 范围 约 1 AU (日 地 平均 距离 ， 
1AU 二 149597870.691km) 采样 154 次 . 采样 间隔 如 
图 10 所 示 多 

从 图 10 可 以 看 出 , 经 度 间 隔 集中 在 4.5°, 浮动 
范围 也 很 小 , 可 以 认为 在 经 度 上 是 等 间隔 采样 的 . 纬 
度 上 间隔 是 非 均 匀 的 , 在 低 纬 上 间隔 大 , 高 纬 上 间隔 
小 , 赤道 上 间隔 最 小 . 从 图 10(c) 所 示 的 半径 间隔 分 
布 可 以 看 出 , 离 太阳 越 远 , 则 半径 间隔 越 远 , 且 呈 线 
性 趋势 增长 , 增长 斜率 约 为 2.86%, 并 且 在 径 向 的 采 
样 次 数 远 多 于 经 向 和 纬 向 上 的 次 数 . 

根据 上 述 分 析 可 以 看 出 , 太阳 风 数 据 在 经 向 和 纬 
向 的 采样 次 数 接近 , 即 经 纬 向 分 辨 率 基 本 同步 , 而 径 
向 采样 次 数 明显 较 大 , 即 半径 分 辨 率 与 经 纬 球面 分 辨 
率 不 同步 . 若 在 经 纬 径 三 个 方向 采用 相同 的 痢 分 次 
数 , 就 会 失去 部 分 径 向 数据 特征 , 造成 一 定 的 数据 失 
真 . 当 径 向 谢 分 层次 比 球面 上 多 时 , 可 以 解决 此 问题 . 

由 于 八 分 体 数 据 维 度 近似 是 2 x 25 x 23, 即 径 
向 采样 密度 是 球面 的 2 倍 , 所 以 可 采用 SDOG-R 格 
网 对 该 数据 进行 痢 分 , 即 在 SDOG 格 网 的 基础 上 继 
续 在 径 向 细 分 3 次 即 可 满足 分 辨 率 要 求 . 

2.2 可视化 显示 

在 日 心 黄道 坐标 系 下 , 以 太阳 风 粒 子 密度 数据 为 

例 , 粒子 密度 单位 为 m 3:, 由 于 密度 值 数据 量 级 相差 
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太 大 , 并 且 每 个 数值 出 现 的 频率 差别 很 大 , 若 用 一 种 
颜色 的 渐变 描述 , 无 法 表达 数据 特征 , 则 需要 较 多 的 
颜色 描述 数据 特性 . 根据 数值 及 其 频率 分 段 制 作 的 颜 
色 模 板 如 图 11 所 示 . 

图 11 所 示 的 太阳 风 密 度 颜色 对 照 显 示 , 颜色 模 
板 共 有 17 个 等 级 , 颜色 等 级 对 应 不 同 的 密度 数值 范 
围 . 根据 图 11 所 示 的 颜色 模板 对 太阳 风 密 度数 据 进 
行 可 视 化 处 理 , 结果 如 图 12 和 图 13 所 示 . 图 12 给 出 
的 是 原始 数据 , 原始 太阳 风 数 据 约 为 2.34GB, 经 数 
据 处 理 提取 出 某 时 刻 的 密度 数据 为 67.496 MByte). 
图 13 给 出 的 是 剖 分 后 的 数据 , 痢 分 层次 LOD 为 3， 
4, 5 时 , 数据 量 分 别 为 0.989 MByte, 4.062 MByte， 
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图 10 采样 点 经 度 (a)、 纬度 (b) 和 半径 (c) 的 间隔 分 布 
Fig.10 Longitude (a), latitude (b) and radius (c) 


interval distribution of sampling data 
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图 11 密度 颜色 对 照 (单位 m3) 


Fig. 11 Comparison drawing for density and color (Unit m3) 


图 12 原始 太阳 风 绝 对 密度 数据 
Fig.12 Original solar wind density data 


11.855 MByte. 可 以 看 出 , 原始 数据 在 中 心 非常 密集 ， 
外 围 稀 玻 , 剖 分 后 的 数据 比较 均 义 , 能 完整 表达 太阳 
风 的 整体 特征 ， 数 据 分 辨 率 随 LOD 值 的 增加 而 增 
大 . 可 视 化 速度 快 , 基本 达到 实时 , 原始 数据 要 加 载 
约 40s 以 上 , 而 放 分 后 LOD=4 数据 只 需要 3s 即 可 
绘制 , 并 且 放大 和 缩小 操作 非常 流畅 . 

如 果 需 要 获取 基于 某 个 特征 或 基于 某 个 区 域 的 
所 有 数据 , 必须 浏览 所 有 文件 才能 提取 对 应 数据 . 在 
信息 巨大 的 情况 下 , 采用 这 种 方式 检索 数据 费时 费力 
且 容 易 遗 漏 , 通过 对 每 个 网 格 进行 编码 , 将 网 格 、 属 
性 数据 、 坐 标 、 时 间 等 信息 统一 组 织 , 可 以 快速 定位 
所 需要 的 信息 . 在 传统 的 检索 中 , 基于 数据 逐 行 查找 ， 
而 基于 该 数据 组 织 方 法 ， 则 将 数据 分 区 分 块 网 格 化 ， 
可 以 层 层 递 进 地 快速 定位 , 只 查找 其 所 在 分 区 , 其 他 
分 区 可 以 不 进行 查找 . 

通过 放 分 和 编码 将 大 基数 据 进 行 有 效 组 织 : 第 
一 , 可 以 提供 多 分 辩 率 层次 数据 , 提高 数据 检索 效率 ; 


第 二 , 可 以 将 空间 位 置 与 属性 数据 关联 起 来 ; 第 三 , 可 
以 有 效 地 减少 元 余数 据 量 , 提高 数据 的 可 用 性 . 


采用 适应 性 SDOG-R 格 网 模型 对 日 地 空间 进行 
立体 剖 分 , 并 提出 了 相应 的 编码 方案 ， 经 试验 验证 ， 
SDOG-R 格 网 模型 不 仅 能 够 很 好 地 解决 日 心 附近 格 
网 过 密 问题 , 还 可 以 满足 径 向 分 辨 率 与 球面 分 辨 率 同 
步 和 不 同步 两 种 需求 ， 同 时 能 够 提供 多 种 分 辨 率 层 
次 的 数据 , 有 效 提高 检索 效率 , 适用 于 不 规则 采样 数 
据 的 痢 分 . 但 是 , 本 文采 用 的 方法 在 半径 向 为 等 间隔 
的 , 如 果 能 够 实现 径 向 非 等 间隔 的 递归 剖 分 , 将 会 更 
加 贴 合 太阳 风 数 据 的 物理 特征 . 

该 模型 不 局 限于 太阳 风 数 据 的 痢 分 , 理论 上 对 类 
似 采 样 的 空间 数据 都 可 适用 . 在 剖 分 的 网 格 中 放置 不 
同 的 数据 , 即 可 实现 多 种 空间 物理 数据 的 组 织 管理 . 
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